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ABSTRACT 

This  paper  contains  an  analysis  of  the  complex  set  of  periodic  solutions 
which  may  occur  in  a  fluid  filled  vessel  of  rectangular  cross  section.  A 
previous  analysis  by  Verma  and  Keller  [2]  -‘found  only  simple  eigenvalues  for 
the  linearized  problem.  It  is  shown  herein  that  at  critical  values  of  the 
vessel  aspect  ratio  double  eigenvalues  occur.  Eight  non-linear  solution 
branches  are  emitted  from  these  double  eigenvalues.  The  solutions  along  the 
various  branches  are  derived,  and  the  results  displayed  graphically.  It  is 
shown  that  irregular  waves  occur  along  some  of  these  branches. 

In  an  interesting  development,  Bauer,  Keller,  and  Reiss  [4],  in  their 
analysis  of  shell  buckling,  showed  that  the  splitting  of  multiple  eigenvalues 
may  lead  to  secondary  bifurcation.  This  theory  is  applied  to  the  non-linear 
standing  wave  problem  herein,  and  it  is  shown  that  secondary  bifurcation  does 
occur  in  the  neighborhood  of  the  double  eigenvalue.  A  perturbation  method  is 
used  to  find  the  secondary  bifurcation  points,  and  the  solutions  along  the 
secondary  branches,  in  the  neighborhood  of  their  respective  branch  points,  are 
found. 

The  neighborhood  around  the  critical  aspect  ratios  is  substantial, 
suggesting  that  secondary  branching  may  occur  in  a  variety  of  vessels  with 
rectangular  cross  section.  This  theory  offers  an  explanation  of  the  irregular 
waves  often  observed  in  the^*sloshing,rof  fluid  in  a  rectangular  vessel. 
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SIGNIFICANCE  AND  EXPLANATION 


This  paper  contains  an  analysis  of  the  complex  set  of  periodic  solutions 
which  may  occur  in  a  fluid  filled  vessel  of  rectangular  cross  section.  A 
previous  analysis  by  Verma  and  Keller  [2]  found  only  simple  eigenvalues  for 
the  linearized  problem.  It  is  shown  herein  that  at  critical  values  of  the 
vessel  aspect  ratio  double  eigenvalues  occur.  Eight  non-linear  solution 
branches  are  emitted  from  these  double  eigenvalues.  The  solutions  along  the 
various  branches  are  derived,  and  the  results  displayed  grapahically.  It  is 
shown  that  irregular  waves  occur  along  some  of  these  branches. 
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ON  THE  SECONDARY  BIFURCATION  OF  THREE 
DIMENSIONAL  STANDING  WAVES 

Thomas  J.  Bridges 

1.  Introduction 

In  a  rectangular  basin  with  vertical  sides,  it  is  well  known  that,  if  the  amplitude 
is  very  small,  the  wave  field  in  the  basin  will  be  a  set  of  elementary  cosine  waves  in  the 
two  horizontal  dimensions  with  a  particular  natural  frequency,  which  is  merely  a  function 
of  the  basin  geometry  and  gravity  (viscosity,  and  surface  tension,  etc.,  being 
neglected).  This  is  the  classical  linear  standing  wave  solution  [1,  Iamb,  p.  364].  As 
the  amplitude  of  the  wave  field  is  increased,  the  number  of  solutions  and  their  type  is 
increased  dramatically. 

Considering  a  cuboidal  container,  the  relevant  parameters  are  the  dimensionless 
depth,  6,  and  the  aspect  ratio,  £,  of  the  horizontal  cross  section.  As  the  amplitude 
of  the  %rave  field  is  increased,  these  two  parameters  govern  the  type  and  multiplicity  of 
solutions  which  arise.  However,  there  is  two  distinct  types  of  secondary  branching  which 
take  place  in  this  problem.  One  is  due  to  the  parameter  £  and  Is  only  weakly  dependent 
on  £  (In  fact,  this  bifurcation  involves  branching  away  from  standing  waves  to  other 
types  of  waves  (sometimes  travelling  waves)),  and  the  other  is  due  to  the  parameter  £ 
and  is  only  weakly  dependent  on  £.  In  both  cases  the  other  parameter  affects  the 
secondary  bifurcation,  and  the  Interaction  can  be  complex,  but  the  two  phenomena  remain 
when  the  other  parameter  is  excluded.  This  enables  the  two  secondary  bifurcations  to  be 
isolated  and  studied  individually,  ttie  bifurcation  which  is  due  to  £  and  remains  when 
the  flow-field  is  two-dimensional  will  be  considered  in  a  separate  study.  The  purpose  of 
this  paper  is  to  study  the  multiple  and  secondary  bifurcations  which  arise  due  to  the 


effect  of  ? ,  the  basin  aspect  ratio.  A  rectangular  vessel  with  vertical  sides  and 
infinite  depth  is  adopted  to  illustrate  this  phenomena. 


The  lowest  natural  frequency,  of  such  a  basin,  and  the  solutions  in  the  neighborhood 

of  this  point  have  been  found  by  Verma  and  Keller  [2] .  These  solutions  were  found  using  a 

perturbation  procedure  in  the  amplitude,  e,  and  they  determined  the  effect  on  the 

natural  frequency,  of  the  amplitude,  to  O(e^) .  However,  by  concentrating  on  a  single 

mode,  this  analysis  missed  the  complex  interactions  which  may  take  place  between  some 

modes  and  the  secondary  bifurcation  which  also  occurs  on  some  primary  branches.  In 

Section  2,  an  analysis,  similar  to  Verma  and  Keller,  is  used  to  determine  the  complete  set 

of  primary  bifurcation  points  and  the  solutions  along  the  simple  primary  branches.  This 

I/- 

analysis  shows  that  each  of  the  primary  bifurcation  points  occurs  at  on  “  X  *  ,  where 

*  u  nif  n 

1  -  ir(m2  +  52n2)  1/2  (1.1) 

ni/  n 

where  m  and  n  are  the  mode  numbers,  respectively,  in  the  x-  and  s-directions.  For 
particular  combinations  of  £ ,  m  and  n,  multiple  eigenvalues,  of  multiplicity  two, 
will  occur.  For  example,  if  ?  “  1  (a  square  basin),  every  pair  (m,n),  m  *  n  is  a 
double  eigenvalue.  If  £  -  V2  »  then  X^4  -  *2,2'  ®tc*  In  Cact'  Cor  5  equal  to  any 
rational  fraction  there  is  a  set  of  double  eigenvalues.  The  physical  reason  for  the 
appearance  of  double  eigenvalues  is  quite  obvious.  In  a  square  vessel,  for  example,  the 
linear  sinusoidal  wave  field  with  2 -wavelengths  in  the  x-dlrection  and  1 -wavelength  in  the 
z-direction  will  share  the  same  natural  frequency  with  a  wave  form  which  has  1-wavelength 
in  the  x-direction  and  2  wavelengths  in  the  z-direction.  Therefore,  although  the 
eigenvalue  is  double,  the  eigenfunctions  are  linearly  independent.  This  is,  in  fact,  a 
ramification  of  spatial  symmetry,  a  subject  which  is  discussed  in  a  more  general  context 
by  Sattinger  [3]. 

In  Section  3,  an  analysis  of  these  double  primary  bifurcation  points  is  performed. 

It  is  found  that  from  each  of  these  points  there  extends  8  branches  (4  of  these  branches 
correspond  to  -e,  and  are  reflections  of  the  +e  branches).  Of  the  four  +e  branches 


-2- 


two  correspond  to  the  two  linearly  independent  eigenfunctions  taken  separately,  and  the 
other  two  correspond  to  mixed  modes  involving  both  of  the  eigenfunctions .  These  mixed 
modes  combine  to  produce  irregular  surface  wave  fields  which  for  some  cases  appear  to  be 
spatially  random  to  the  uninitiated  observer.  Figures  depicting  solutions  of  this  type 
are  given  in  Section  3. 

In  an  intersting  discovery,  Bauer,  Keller,  and  Reiss  [4]  observed  that  as  a  multiple 
bifurcation  point  is  "split”  (by  varying  an  auxiliary  parameter  (in  this  case  £))  into 
primary  bifurcation  points  a  secondary  bifurcation  may  occur.  At  the  double  eigenvalue, 
which  arises  in  the  standing  wave  problem,  there  are  four  +e  branches.  Upon  splitting 
this  point,  two  of  these  branches  become  the  simple  primary  branches.  Due  to  continuous 
dependence,  the  other  two  solutions  cannot  simply  vanish  when  £  is  perturbed  away  from 
£q,  the  location  of  the  double  eigenvalue,  and  therefore  they  slowly  depart  by  creeping 
up  the  split  primary  branches  in  a  manner  which  differs  from  problem  to  problem.  Bauer, 
Keller,  and  Reiss  developed  a  perturbation  method  which  enabled  them  to  analyze  this 
phenomena  in  the  buckling  of  thin  shells.  Subsequently,  this  theory  has  been  successfully 
applied  by  Mahar  and  Matkowsky  (5]  to  a  model  biochemical  reaction  problem,  Matkowsky, 
Putnick,  and  Reiss  [6]  to  the  buckling  of  rectangular  plates,  and  has  been  extended  to  the 
bifurcation  from  triple  eigenvalues,  which  results  in  secondary  and  tertiary  bifurcation, 
by  Reiss  [7] . 

In  Section  4  the  theory  of  Bauer,  Keller,  and  Reiss  is  used  to  show  that  secondary 
bifurcation  occurs  in  the  neighborhood  of  the  double  points  which  arise  in  the  problem  of 
three-dimensional  standing  waves.  The  secondary  bifurcation  points  are  found  and  the 
solutions  along  them  are  derived.  It  is  found  that  mode  Jumping  occurs  in  the  shift  from 
the  primary  to  the  secondary  branches  and  the  complex  interaction  of  the  two  modes  in  the 
neighborhood  of  the  secondary  bifurcation  point  produces  irregular  wave  forms. 

These  mixed-mode  solutions  offer  an  explanation  of  some  of  the  irregularities  that 
are  often  observed  in  experiments.  For  example,  in  the  experiments  of  Taylor  [8],  on  the 
highest  two  dimensional  standing  waves,  irregular  three-dimensional  "crown  waves"  were 


observed  forming  on  top  of  the  two-dimensional  standing  waves.  His  experimental  basin  had 
an  aspect  ratio  of  C  ~  2.0,  and  the  crown  waves  look  like  mixed-mode  solutions  discussed 


herein . 


*  wv 


2.  Primary  Bifurcation  from  Single  Eigenvalues 


The  fluid  in  the  basin  is  taken  to  be  inviscid,  irrotational ,  and  absent  of  surface 


tension-  The  problem  can  then  be  expressed  in  terms  of  a  potential  function  and  the  wave 


height-  Length  scales  are  rendered  dimensionless  by  2a,  the  tank  length  (see  Figure  1), 


with  the  exception  of  the  z-direction  where  2b,  the  tank  width  is  used,  and  g,  the 


acceleration  of  gravity,  and  the  frequency  arc  used  for  the  time  scales.  The  governing 


equations  and  boundary  conditions  are 


2  *  2  *  2  * 

♦  ^-f-  ♦  52  -  o 

3x^  3y  1zl 


for  -  V»  <  «  <  Vo  »  -  4  < 


•  <  y  <  n  (x,z,t) 


-  0  at  x  -  *  V2 


(2.2a) 


0  at  z  ■  ±  V2 


(2.2b) 


•*  0  as  y  ♦  - 
dy 


(2.2c) 


“  lt“  +  2  t(*x)2  +  ( *y)2  +  e2(*z,2]  +  T1# 


•0  on  y  “  n 


u  3n  .*  3n  ,2  *  3n  3$  * 

sT  +  *x  IT  +  *  *z  IT  ~  sT  “  0  on  y  “  n 


$  *=  n  -  0  is  a  solution  for  all  values  of  £.  This  is  the  basic  state.  Periodic 


solutions,  of  period  u  ,  which  bifurcate  from  this  state  are  sought.  The  potential 


function  and  the  wave  height  are  assumed  to  have  an  expansion  of  the  following  form. 


4>  (x,y,z,ti  e)  «  e*1  (x,y,z,ti  5)  +  e  $  (x,y,z,t>  £)  +  .  •  • 


(2.5a) 


‘  1  .**  fc'*  ,*•  *•  's  “«  ' 


-  •  •.'-y  .  v-v  ■  • 


tv  * 


v.-v-v-  >  *  vN 

,  v-v-.-.v.y 

.•v.-Vsl  -.- -A  o  .  \\?l  v\\.A 


n  ( x, z , t ;  5 ; e )  «  er^  (x,z,t; t)  +  £  h2(x,z,ti£)  +  .  .  . 


(2.5b) 


ID  (?>£)  =  aQ ( e>  +  £0^5)  +  £  e>2 C C)  +  .  .  . 


(2.5c) 


and  e,  a  measure  of  the  wave  amplitude,  is  chosen  to  be  proportional  to  the  first  term 
in  (2.5b).  Before  substitution  of  (2.5)  into  the  boundary  value  problem,  the  free  surface 
boundary  conditions  (2.3)  and  (2.4)  are  expanded  in  a  Taylor  series  about  y  »  0, 


F(x,y,z,t ;£>£) 


3f i  .*  .  1  32F.  *2 


n  *Ty\\  +  2“2  *  + 

y»0  1 y*0  3y  1 y*0 


Substitution  of  (2.5)  into  (2.1)  -  (2.4),  the  use  of  (2.6),  and  the  setting  to  zero  of  the 
terms  proportional  to  each  power  of  e  results  in  the  following  sequence  of  boundary 
value  problems, 

sY  »Y  ,  aY 

- 1  + - 1  +  ^ - 2-  -  o 

3x2  3y2  3z2 

in  -y<x<j,  -»<y<0  (2.7) 


0  at  x  »  t  V2 


(2.(?a) 


*  0  at  *  -  ±  V2 


(2.8b) 


♦  0  as  y  ♦  - 


o.  -xr-2-  +  n.  «  R.  at  y  -  0 


3y 

34 

30  It  * 


3n  34 

j.  -r— 2-  -  ■» — —  »  S  .  at  y  -  0 
0  at  3y  j  1 


(2.8c) 


(2. 10) 


for  j=1,  2,  ...<  R1  -  S1  «  0,  and  Rj  and  Sj  for  j  >  1  are 


functions  of  terms  of 


previous  order  only,  and  are  given  in  Appendix  I.  Since  (2.7)  -  (2.10)  for  j  >  1  are 
inhomogeneous  eigenvalue  problems,  they  are  soluable  if  and  only  if  they  satisfy  the 
Fredholm  alternative.  The  alternative  necessary  for  this  and  the  later  developments  is 


given  in  Appendix  II. 


Wte  solution  of  (2.7)  -  (2.10)  for  j  ■  1  is 


“  cos  c^x  cos  Bn*  sin  t 
- 1  ny 

•  o.  e  '  cos  a  x  cos  8  *  cos  t 
10  m  n 

where  a  “  8  ”  nr,  x  •  x  ♦  V>  ,  z  «  z  *  V>  ,  and 

m  n  * 

o  2  -  X  -  (a  2  +  C28  2)  1/2 
0  m,n  m  n 

At  the  next  order,  the  solvability  condition  requires  that  c^  •  0,  then 

n.  “  (B  +  C  cos  2t)  cos  2a  x” 

2  21  21  01 

*  (B__  +  C,_  cos  2t)  cos  28  z” 

22  22  n 

+  (B_,  +  C„,  cos  2t)  cos  2a  x  cos  28  z 

23  23  m  n 


2amy  -  258ny 

*2  -  [A20  +  A21e  cos  2amx  +  A22e  cos  28n*]  sin  2t 


(2.11) 


(2.12) 


(2.13) 


(2.14) 


A20  “  *  V8 


21  2 
4a  (a  -  2a  ) 
0  m  0 


“22  4a0(C8n  -  2a02) 


21  _  2 
8°0 


(2.15) 


(2.16a) 


(2.16b) 


(2.16c) 


(2. 17a) 


(2.17b) 


0  /8 


(2.17c) 


23 


C21  “  "  2O0A21 


<0  4  +  C28  2 ) 
0  n 


(2. 17d) 


8  0. 


-  -  2 0. A 
22  0  22 


<O04  +  am2) 


Bo, 


C23  “  -  °0  /8 


(2. 17e) 

(2. 17f ) 


After  some  manipulation,  the  solvability  condition  at  the  next  order  results  in  an 


expression  for  the  amplitude  correction  to  the  natural  frequency. 


5  ,,4  4.4  ,4.4 

a0  3  °m  +  5  6n  >  5  6n 

°2 - 8 - ~~~3 -  +  .  “  2 


(2.18) 


32  a. 


4o0(2o0  -  «B)  4d0(2o0  -  6Bn) 


By  choosing  a  »  8  »  1  and  5  ”  1/L,  this  result  is  in  agreement  with  that  found  by 

m  n 


Verma  and  Keller  for  the  first  mode  when  the  depth  is  infinite.  The  expression  (2.18)  is 
negative  definite  and  this  can  be  shown  by  recasting  (2.18)  as 


10  8  4  2  4 

-4o.  +  258 o.  +  2a  a,  +  356  a 

r  0 _ n  0 _ m  0 _  n  m  i 

°2  “  °0  A  7  J 


32a04  (2o02  -  5Bn) 


—4c.  +  2a  a  8  +  2?4B  40  2  +  3a  £4B  4 

_  r  0 _ mO _ n  0 _ m  n  1 

°n  L  a  7  \ 


32a  4(2o  2  -  a  ) 
0  0m 


(2.19) 


with  no  loss  of  generality,  choose  56  «  ra  where  r  e  [0,«].  Then  (2.19)  can  be 

n  m 


expressed  as 


For  subcritical  bifurcation  it  is  required  that 


2/1  +  r  (1 


2(1  +  r2)2)  +  r (3  +  2(1  +  r2)2)  <  0 


(2.20a) 


and 

2/1  +  r2(r4  -  2(1  +  r2)2)  +  (3r4  +  2(1  +  r2)2)  <  0  (2.20b) 

for  any  r  e  [0,°°].  After  some  algebraic  manipulation  (2.20a)  can  be  recast  as 

-  [4  +  3r2  +  40r4  +  92r6  +  64r®  +  12r10]  <  0  (2.21a) 

and  (2.20b)  can  be  recast  as 

-  [2  +  12p  +  8p2  +  4p3  +  5p4  +  2p5]  <  0  (2.21b) 

where  p  =  J\  +  r2  -  1,  which  satisfies  p  e  [0,«)  when  r  e  [0,-).  The  pair  (2.21a,b) 
proves  the  conjecture:  the  bifurcation  from  all  the  simple  eigenvalues  is  subcritical. 

The  natural  frequency  of  the  standing  wave  decreases  as  the  amplitude  increases.  There  is 

also  the  interesting  feature  that  | o2 |  increases  with  increasing  m  and  n.  Therefore 

* 

the  magnitude  of  the  slope  in  the  ui  -e  plane  decreases  with  increasing  mode  number 

resulting  in  an  intersection  of  the  higher  mode  branches  with  the  lower  mode  branches. 

Figure  2  is  a  plot  of  the  natural  frequency  as  a  function  of  amplitude  for  5  -  .25 

and  m,  n  ranging  over  1/  2,  and  3.  The  dashed  lines  correspond  to  modes  with  equal 

indices  (m  -  n).  This  is  an  almost  two  dimensional  field  and  the  9  modes  are  coalescing 

around  the  3  two-dimensional  values.  In  Figure  3,  for  the  same  range  of  m  and  n,  and 

5  =  1.1,  the  modes  are  plotted.  Here  they  are  more  spread  out  with  the  mode  intersections 

being  quite  obvious.  In  Figure  4,  a  typical  non-linear  wave  on  a  branch  bifurcating  from  a 

simple  eigenvalue,  is  plotted.  The  vessel  is  square  in  cross-section  (£  »  1),  e  “  0.15, 

and  m  =  n  =  1.  This  is  often  referred  to  as  the  "slosh"  mode.  Figure  5  shows  a  time 

series  for  the  wave  height  in  the  left  front  corner  (x  ■  z  «  0),  of  the  tank  depicted  in 

t  3 

Figure  4.  Figure  4  is  a  "snapshot"  corresponding  to  —  -  —  ±  k,  k  -  0,1,2,  .... 


3.  Primary  Bifurcation  From  the  Double  Eigenvalues 

At  critical  values  of  the  aspect  ratio  of  the  vessel  cross-section,  a  set  of  double 
eigenvalues  will  occur.  When  £  «  1  there  is  a  double  eigenvalue  for  any  pair  (m,n) 
such  that  m  *  n,  and  when  £  =  V2  (m,n)  =  (1,4)  and  (m,n)  =  (2,2)  is  a  double 
eigenvalue,  ad  infinitum.  For  €  equal  to  any  rational  fraction  there  is  a  set  of  double 
eigenvalues.  It  is  expected  that  the  nature  of  the  solutions  will  be  similar  for  the 
various  critical  values  of  £.  To  elaborate  this  phenomena,  an  analysis  of  the  solutions 
emitted  from  the  set  of  double  eigenvalues  for  a  square  vessel,  5  »  1,  is  performed.  In 
this  case  the  bifurcation  points  are  given  by 

°0  =  (“m2  +  6n2)  ^  (3‘1) 


The  potential  function  and  wave  height  are  again  expanded  in  a  perturbation  series 
in  e , 

♦*(x,y,z,tjC >e)  =  (x,y,z,t»S)  +  e2$2(x,y,z,t>£)  +  .  .  .  (3.2a) 

n  (x,z,t;£,e)  =  (x,z,t;C)  +  e2n2(x,z,t»E)  +  .  .  .  (3.2b) 

w*(C,e)  =  o0(5)  +  e2^^)  +  e2o2(C)  +  .  .  .  (3.3) 

However,  the  symbols  here  should  not  be  confused  with  those  in  Section  2.  Outside  of  the 
obvious  similarities,  the  results  in  Section  2  are  distinct  from  the  results  contained 
herein.  In  fact  they  are  the  complement  of  each  other. 

Substituting  (3.2)  and  (3.3)  into  the  governing  equations  and  boundary  conditions 
(2.1)  -  (2.4)  and  equating  terms  proportional  to  like  powers  of  e  to  zero  results  in  the 
sequence  of  boundary  value  problems  given  in  (2.7)  -  (2.10)  but  with  the  R  ^  and  S  ^ 
differing  upon  substitution. 

with  R1  *=  S1  »  0,  and  the  fact  that  there  is  a  double  eigenvalue  with  a  two- 
dimensional  null  space,  the  first  order  solution  is 

n.  =  (A. .cos  a  x  cos  8  z  +  A, .cos  8  x  cos  a  z]  sin  t  (3.4) 

1  iim  n  n  m 


and 
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♦«  •  a.  [A,, coo  a  x  cos  0  z  ♦  A. .cos  0  x  cos  a  z]  cos  t  (3*5) 

1011  m  n  12  n  m 

where  a,,  is  givsn  by  (3.1)  and  e  is  chossn  to  be  proportional  to  the  first  order  wave 

height.  This  rssults  in  a  relationship  betwesn  and  A12» 

A,,2  ♦  A122  -  1  (3.6) 

The  relative  magnitudes  of  Ajj  and  A12  are,  however,  undertermined  at  this  order. 

The  problems  for  j  >  1  are  inhomogeneous  eigenvalue  problems  and  must  satisfy  the 

solvability  condition  given  in  Appendix  II.  As  the  nullspace  at  the  eigenvalue  is  of 

dimension  two,  there  will  be  two  solvability  conditions  at  each  order.  Substituting 

♦  .  ,h .  ,etc.  into  Rj  and  S2,  and  applying  the  solvability  condition  results  in 

Oj  ■  0,  and 

n2(x,s,t)  -  (B21  +  C21cos  2t)  cos  2amx 

+  (BJ2  C22cos  2t)  cos  2a#S 

+  <B23  +  C23°°*  2t*  co“ 

+  (B,,  +  C-.cos  2t)  cos  26  i 
24  24  Ti 

+  (B  +  C,ccos  2t)  cos  (a  +  6  )x  cos  (a  ♦  6  )x 
2d  23  n  n  n  n 

+  (B  +  C  eos  2t )  [costa  ■¥  6  )x  cos ( a  -  6  )*  + 

20  20  m  n  n  n 

cost  a  -  8  )x  cos(o  ♦  6  )i] 
n  n  n  n 

+  (B  ♦  C__  cos  2t)  costa  -  6  )x  costa  -  6  )* 

2/2/  a  n  nn 

+  IB,.  +  C„  cos  2t)  cos  2a  x  cos  26  z 
2o  20  m  n 

+  (B  +  C  cos  2t)  cos  28  x  cos  2a  x  (3.7) 

2i»  2^  n  n 


♦2(x,y,z,t)  -  [a20  +  A21cos  2aBx  exp  Uc^y) 
+  AJ2  cos  20^2  exp  (2amy) 


+  A,,co8  26_x  exp  (26_y) 


+  A.. cos  20  z  exp  (20  y) 
n  n 

+  A  cos  (a  +  8  )x  cos(a  +  8  )*  exp  [/2  (a  ♦  8  )y] 
43  m  n  on  n  n 


A_,[cos  (a  8  )x  costa  -  6  )z  + 
2o  in  n  n  n 


+  cos  (a_  -  8  )x  cos  ta  +6  )z]exp  [/2a.  yj 
m  n  m  n  u 

+  A  cost  a  -  8  )x  cos  (a  -  6  )z  exp  [/2|a  -  8  |y] ]sin  2t  (3.8) 


The  coefficients  A2 j ,  Bj j <  end  C2 j  are  given  in  Appendix  III. 

The  most  obvious  effect  of  the  double  eigenvalue  is  the  explosion  of  higher  harmonics 
at  the  second  order.  Ibis  is  to  be  contrasted  with  the  solutions  at  simple  eigenvalues 
where  simple  harmonics  appear  at  higher  order  and  only  weakly  affect  the  linear  solution. 
In  the  case  of  the  double  eigenvalue  the  mixed  mode  solutions  also  result  in  an  increasing 


complexity  in  the  spatial  structure,  on  a  finer  and  finer  scale,  as  the  amplitude  is 
increased. 


Substituting  the  known  terms  into  Rj  and  S3  and  applying  the  two  solvability 
conditions  results  in  the  bifurcation  equations, 

[Vl12  +  Vl22  +  2°2]  *11  “  0  (3 


lVl12  +  «1*122  +  202]  *12  “  0 


(3.10) 


which,  along  with  the  normalisation 


*112  +  *122  ’  1 


(3.11) 


form  a  set  of  three  equations  for  the  three  unknowns  A^,  A12,  and  a^.  The 
coefficients  a^  and  a2  are  given  by 


5  4  4  4 

.  q0  i  3<am  +  6n  > _ Bn 

4  16o03  2O0(2O02  -  Bn)  2O0(2°02  -  V 


( 3 • 12) 


.  *.*  *  ■. 


•  •  .  .  \  •  V  .  „■ 

V.--V  • 

VV  A 


y  \vy.\y.--Vv'  v-\*-y-' 


2 


^10  -  3/2 i  °0 
4-/2 


3(a  +  0  *) 

\  u  m  n 

J  4 


(am  -  V 


160. 


4o0(4o0 


/2(a  +  e„)) 

m  n 


(a  +  6  ) 

m  n 


4o0(4ofl2  -  /2|an  -  Sj) 


The  three  equations  (3*9)  -  (3.11)  has  the  following  set  of  solutions 
Case  It  Aq2  ”  0,  AQ1  «  ±1#  o2  -  -a.j/2 

Case  II:  AQ1  -  0,  A02  -  ±1,  Oj  -  -aj/2 

Case  III:  A01  -  i  —  ,  A  -  i—  ,  0*  -  -(a1  +  a2>/4 

/2  /2 


(3.13) 


Case  I  and  II  are  pure  modes  corresponding  to  the  modes  of  the  simple  eigenvalues  that 
coalesce  to  form  the  double  point.  They  share  the  same  natural  frequency  and  are  spatially 
symmetric.  When  the  relevant  parameters  are  substituted  the  amplitude  correction  to  the 
natural  frequency,  o2  «  “  'J  a1'  a9re*B  with  the  correction  found  for  the  simple 
eigenvalues  (Eqn.  (2.18)).  In  Section  2  it  was  found  that  this  expression  is  negative 
definite,  therefore  a1  >  0  for  all  positive  integers  m  and  n.  Case  III  involves  mixed 

modes.  The  leading  terms  are  proportional  to  the  sum  and  difference  of  the  two 

M 

eigenfunctions.  The  amplitude  correction  to  the  mixed  mode  solutions,  ,  differs  from 

the  pure  mode  case  by  an  amount  ( a  ^  -  a2 )/4.  It  has  not  been  proved  that  this  expression 
is  negative  for  all  parameter  values,  but  it  appears,  from  numerical  evaluation,  that  this 
is  in  fact  the  case.  It  is  therefore  conjectured  that  the  bifurcation  of  the  mixed  modes 
is  always  subcrltical. 

Figures  6  and  7  display  the  bifurcation  from  double  eigenvalues  for  a  vessel  with  a 
square  cross-section.  In  Figure  6  the  dashed  lines  correspond  to  the  simple  eigenvalues 
I,  1  and  X2  2,  and  the  solid  lines  are  emitted  by  the  double  point  X^  2*  There  are 
two  pure  modes  corresponding  to  the  upper  branch  (labelled  p)  and  two  mixed  modes  which 


Figure  6.  Bifurcation  diagram  of  simple  and  double  eigenvalues 

for  £  =  1.0  and  m  and  n  ranging  over  1  and  2. 

The  dashed  lines  correspond  to  the  A  and  A. 

1 r 1  2t2 

branches.  The  upper  branch  emitted  by  the  double  eigen- 
supports  two  pure  modes  (labelled  p) ,  and  the  lower  branch 
supports  two  raided  modes  (labelled  m) . 


Bifurcation  of  simple  and  double  eigenvalues  for  £  *  1.0 
and  (m,n)  ranging  over  1  and  3.  The  dashed  lines  cor¬ 
respond  to  the  A  and  A  branches.  The  upper 
1  f  X  -i  9  J 

branch  emitted  by  the  double  eigenvalue  supports  two  pure 
modes  (labelled  p) ,  and  the  lower  branch  supports  two  mixed 
modes  (labelled  m) . 


share  the  natural  frequency  on  the  lower  branch  (labelled  m) .  Figure  7  shows  a  similar 

result  with  m  and  n  ranging  over  1  and  3.  The  dashed  lines  correspond  to  X f  j  and 

X,  ,,  and  the  solid  lines  correspond  to  the  four  branches  emitted  by  the  double  point. 

In  Figures  8,  9,  and  10  the  wave  fields  corresponding  to  Figure  7  are  plotted.  In 

Figure  8,  for  m  •  1,  n  *■  3,  and  e  ■  .1,  the  pure  mode,  with  A^  =  1.0  and  A-|2  “  0 

is  plotted.  This  corresponds  to  the  upper  branch  emitted  from  the  double  eigenvalue  in 

Figure  7.  In  Figures  9  and  10,  the  results  along  the  lower  (mixed  mode)  branch  emitted 

from  the  double  eigenvalue  are  plotted  for  e  «  .1.  In  Figure  9  A^,  •  —  and 

1  ^ 

A, _  *  —  ,  whereas  in  Figure  10  A  ■  —  and  A,  ■  -  —  .  The  mixing  of  the  modes 

12  a  11  /s  12  a 

produces  quite  irregular  solutions.  Figures  11a,  b,  c  show  the  time  series  for  the  left 
front  corner  (x  =  z  =  0)  of  the  vessels  depicted  in  Figures  8,  9,  and  10.  The  snapshots 
in  Figures  8,  9,  and  10  correspond  to  the  times  ^  “  |  ±  k,  k  «  0,1,2,  ....  Figure  11a 
corresponds  to  the  pure  mode  in  Figure  8  (A^  •  1,  A12  »  0).  Figure  11b  corresponds  to 
the  mixed  mode  in  Figure  9  (A^  «  1//2,  A^2  “  1//2).  Figure  11c  corresponds  to  the  mixed 

mode  in  Figure  10  (A^  «  1//2,  A12  -  -1//2).  Reference  to  eqn.  (3.4)  shows  that  n1  “  0 

-  -  2 
when  x  =  z  =  0  and  a^ 1  -  “A12-  Hence  the  linear  effects  in  Figure  11c  are  of  0(e) 

and  proportional  to  cos  2t,  with  the  nonlinear  effects  ccming  in  at  0(e3).  Although  the 

results  in  Figures  11a,  b,  and  c  differ  markedly,  all  three  correspond  to  the  same  set 

of  basic  parameters. 

Examples  for  other  mixed  mode  solutions  in  a  square  vessel  are  given  in  Figures  12  and 
13.  Figure  12  is  of  a  mixed  mode  with  m  “  2,  n  «  4,  and  Figure  13  is  of  a  mixed  mode 
with  m  =  1,  n  »  5. 
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Figure  8. 


,  ~e  solutions  along  the  p 

Spatial  variation  o°  t  k  -  0,  1,  2,  .. 

branch  in  Figure  7  at  t/2ir  /  0.1, 

The  parameters  are  £  -  1.,  #>  -  1<  n 

«  _ J 


All“  l" 


and 


A  —  0. 

*12 


Spatial  variation  of  one  of  Re  solutions  along Rhe^m 

r^-riri  i  -  a.  e = °-1- 

a  =  1/^,  and  A  ,  =  V«£* 


Figure  9. 


Figure  10.  Spatial  variation  of  the  other  solution  along  the  m 

branch  in  Figure  ^  at  t/2t  -  3/4  i  k,  k  =  0,  1,  2,  .... 
The  parameters  are  £  =  1,  m  =  1,  n  *  3,  e  *  0.1, 

R11  =  304  A12  “ 


Figure  11a. 


Time  series  of  the  wave  height  in  the  left  front  comer 
(x  =  z  =  0)  of  the  vessel  in  Figure  8.  The  parameters 
are  £  =  1.,  m  =  1,  0*3,  £  “  0.1,  A^  =  1. ,  and 


-23- 


EiSg 


branch  esnitted 

.e  o£  the  wave  Com  orw  mixed  ^  ^  «  i,  e  -  0.1. 

-rbl:rra:f:.^-  vr“-- 

“a  *u  >2  »  p:od»c..  .  **•«“  a««*  <*  s»” 

is  a  multiple  of  m  P 
in  the  mixed  mode. 


Held  illuetr.tind  the 


The  Secondary  Bifurcation  Points 


In  this  section,  solutions  which  bifurcate  from  the  primary  branches  found  in  Section 
2  are  sought.  Therefore  a  perturbation  is  added  to  the  knowr.  primary  branch  solution. 


=  e  4  +  4> 


(4.1a) 


e  h  +  n 


(4.1b) 


The  expressions  (4.1!  are  substituted  into  the  governing  set  of  equations  and  boundary 
conditions  (2.1)  -  (2.4)  and  linearized.  The  equations  independent  of  $  and  n  were 
analyzed  in  Section  2.  The  linear  problem  for  ♦  and  n  is 

iff  +  iff  +  52  iff  =  0  (4.2 

3x2  3y2  3z2 


0  at  x  =  *  V2 


(4.3a) 


3$ 

87  +  0  as  y 


(4.3b) 


£♦0  at  z  =  *  V2 


(4.3c) 


and  on  y  -  c  h(x,z,t). 


*  in  +  *  ii  in  +  *  ii>  ii 

3t  3x  3x  3x  3x 


*2  f34  3n  3h  36  ,  3* 

+  e  5  1  te  te  +  8?  3z  1  _  3y  =  ° 


*  3*  *  r  34  34  34  34  2  34  34  i  A 

w  t31  +  e  hr-  -r*  +  -r—  -r*  +  £  t-  +  n  *  0 

3t  ^-Sx  3x  3y  3y  s  3z  3z  J 


This  is  a  linear  differential  eigenvlue  problem  with  known  nonconstant  coefficients. 


However,  the  specific  value  of  €  where  the  secondary  bifurcation  takes  place  is 
*  * 

sought.  Therefore  e  is  the  eigenvalue.  Since  e  appears  non-linear ly,  it  is  a  non¬ 
linear,  in  the  parameter,  eigenvalue  problem.  This  does  not  pose  serious  complications  as 
methods  for  the  solution  of  this  type  of  problem  are  known  [9J .  However,  there  is  the 


1  ,  *  *  *  •  '  t  "  •  '  •  *  »  *  •  *  •  ■  »  *  ■ 

•  ■  «  *  •  %  •  m  •  '  ■ 


NV-V 


*  \  .Vv'.V 


further  complication  that  e*  is  responsible  for  the  size  of  the  domain.  Hie  qualitative 
shape  of  the  domain  is  known  since  h(x,z,t)  is  a  known  function,  but  the  precise  amount 

of  h(x,z,t)  is  the  unknown  eigenvalue.  This  is  to  be  contrasted  with  the  original 

•  • 

eigenvalue  problem  in  Section  2  where  u>  was  the  eigenvalue,  e  was  a  variable 
parameter,  and  the  shape  of  the  free  surface,  and  hence  the  domain,  was  an  unknown 
function. 

To  solve  this  eigenvalue  problem,  the  conjecture  of  Bauer,  Keller,  and  Reiss,  that  the 
secondary  bifurcation  occurs  in  the  neighborhood  of  multiple  eigenvalues,  is  used. 

For  brevity,  the  analysis  is  undertaken  in  the  neighborhood  of  5=1  (a  square 
cross-section).  It  is  expected  that  a  similar  analysis  will  hold  in  the  neighborhood  of 
other  values  of  5  at  which  double  eigenvalues  occur. 

It  was  shown  in  Section  3  that  for  5  «  1  there  is  a  double  eigenvalue  for  every  pair 


(m,n)  such  that  m  *  n.  At  the  double  eigenvalue  the  bifurcation  point  corresponds  to 


if  (m2  ♦  n2)  V* 


for  any  (m,n)  such  that  m  *  n,  and  in  the  neighborhood  of  ?  ”  1  this  double  eigenvalue 
splits  into  two  primary  branches 

\  -  v(m2  +  52n2)  1/2  (4.7) 

m,n 
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with  (4.9)  and  (4.10)  the  solutions  on  the  primary  branch,  derived  m  Section  2,  can  be  re¬ 


expressed  in  terms  of 


h  =  hQ  +  yh1  + 


It)  -  it)„  +  ii  u>„  + 
0  2 


(4.11a) 


(4.11b) 


(4.12) 


The  analysis  is  undertaken  separately  on  the  two  branches,  X  and  X  after 

m#n  n ,  m 

splitting.  The  necessary  details  for  the  analysis  along  the  branch  X  is  given,  and 

01/  n 

the  result  only  will  be  stated  for  the  X  branch.  Along  the  primary  branch 

ii/  in 

corresponding  to  X  , 
ni/  n 


X  -  x„  [i  +  — —  u2  +  .  •  .  ] 

m,n  0  L  2  J 

A0 


(4.13) 


where  XQ  is  given  in  equation  (4.6).  The  terms  in  the  expansion  (4.12)  are  given  by 


"o  =  xo  /2 


(4.14a) 


n  ^  .  2  r  0 

“2=7T+bo“o  l— 

4oj 
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4  4  4 

-o  3«\n+6n> 


+  *„2  [' 


4V2“02  -  V  4“o(2“o2  ’  V 


and  the  terms  in  the  wave  height  expansion  (4.11b)  are  given  by 

h„  =  cos  a  x  cos  0  z  sin  t 
0  m  n 

h1  “  <D11  +  °12  008  2t)  008  2amx  +  <D13  +  D14  008  2t)  cos  2®nz 
+  (D,_  +  D  cos  2t)  cos  2a  x  cos  26  z 


(4.14b) 


(4.15a) 


(4.15b) 


f  f  . 


A  .  *v 


where  a  »  m»,  8  “  nx,  x  «  x  +  V2  »  *  ■  z  +  V2  «  an<3 

m  n 


11  „  2  * 

so) 

o 


-b.(w.4  ♦  8  2)  b.B  2 
0  0  n  On 


2(o  -  2u.  ) 

m  0 


13  _  2 

% 


4  2  2 

-b„(u)„  +  a  *)  b  o 

0  0  m  0  m 


2(6  -  2u>  * ) 

n  0 


(4. 16a) 


(4. 16b) 


(4. 16c) 


(4. 16d) 


D1S  '  Vo'8 


(4. 16e) 


D  «  -  b  a)  ‘/8 
16  0  0  ' 


( 4 . 1 6f ) 


end  a  similar  expression  is  found  for  4  in  (4.11a).  Substituting  the  expressions  (4.9), 
(4.10),  (4.11)  into  the  eigenvalue  problem  (4.2)  -  (4.5)  and  postulating  that 


♦  =*  u*,  +  U  $2  + 


n  -  un,  +  u  n2  + 


(4.17a) 

(4.17b) 


expanding  the  free  surface  boundary  conditions  in  a  Taylor  series,  and  equating  terms 
proportional  to  like  powers  of  u  to  zero,  results  in  the  sequence  of  problems. 
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(4.19b) 
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is. ii  i 


for  j  *=  1,  2,  . ... 

In  addition  to  the  obvious  differences  between  the  Rj  and  Sj  appearing  here  and 


those  in  the  previous  sets,  there  is  also  the  addition  of  the  inhomogeneous  term  in  the 


a2*. 


governing  equation  (4.18).  However  “0,  T2  •  and  T3  ”  ”T  — 5” 

3z 


Noting  that  “  0,  and  the  fact  that  uQ  corresponds  to  a  double 

eigenvalue,  the  first  order  solution  is 


n.(x,z,t)  «  [A,, cos  a  x  cos  B  z  +  A. .cos  B  x  cos  a  z] 
1  n  m  n  12  n  ra 


(4.22) 


and 


-  -  -  -  “0  y 

(x,y ,z,t)  ”  ui_  (A.. cos  a  x  cos  B  z  +  A  cos  B  x  cos  a  z]e  cos  t  (4.23) 

1  0  11  m  n  12  n  m 


For  definitness  in  the  normalization r  the  following  relation  is  taken  between  Ajj 
and  A^2» 


An2  +  *i22  '  1 


(4.24) 


with  the  relative  values  of  An  and  A12  to  be  determined  at  higher  order. 

Before  proceeding  to  higher  order,  it  is  noted  that  the  homogeneous  problem  has  a 
double  eigenvalue.  Therefore  two  solvability  conditions  are  needed  at  each  order.  The 
general  form  of  the  solvability  condition,  given  in  Appendix  II,  is  used.  Substitution  of 
the  known  solutions  into  R2  and  S2  the  second  order  solution  is  found  to  be, 


n2(x,z,t)  *  C2oboAi icos  ^ 


♦  [®21  +  C21C°8  2tI  Vl1C°8  2V 


+  [B22  +  C22cos  2t]  b0Ancos  2Bnz 


+  +  <^c°8  2t5  bnA,icos  2amx  008  26nz 


I 


I 


i . 

» 

r- 


+  [B„.  ♦  C_.cos  2t]  b.A._cos  (a  -  8  )x  cos  (a  -  8  )z 

24  24  0  12  n  n  zb  n 


*  [B,.  +  C .cos  2t]  b.A  cos  (a  +  0  )x  cos  (a  +  0  )z 
2D  23  u  12  m  n  m  n 


+  c«bnAiioos  21  ^cos  <a_  “  6_)x  co.  (a  +  0  )z 
2o  u  12  an  ib  n 


+  cos  (a  +  0  )x  cos  (a  -  0  )*] 
in  n  in  n 


(4.25) 


and 


VX'y'*'t>  “  A00b0A11t  +  A20b0A1l“in  21 

2a  y  _  20  y 

+  b„A«.sin  2t  [A_.e  cos  2a  x  +  A,.e  cos  20  z 
011  21  IB  22  Z1 

u02y 

+  A,,e  cos  2a  x  cos  20  z] 

23  m  n 

+  b  A  sin  2t[A  cos  (a  -  0  )x  cos  (a  -  0  )z  exp[/5|a  -  0  | y ] 

0  12  L  24  m  n  m  n  1  n  n' 

+  A. _cos  (a  +  0  )x  cos  (a  +  0  )z  exp  [/2  (a  +  8  )y] 

23  m  n  m  n  an 

+  A  cos  (a  -  8  )x  cos  (a  ♦  0  )z  exp  [/2  un2y] 

2o  an  m  n  u 

+  A  cos ( a  +  0  )x  cos  (a  -  0  )z  exp  w  2y] ]  (4.26) 

26  m  n  m  n  u 

The  coefficients  A2j,  B2j,  and  C2j  are  given  in  Appendix  IV.  Hie  term  C20b0A1  ico®  2^ 
in  the  expression  for  n2  in  (4.25)  is  indicative  of  non-conservation  of  mass. 

Conservation  of  mass  requires  that 

V2  V2 

/  /  n  (x,z,t)dxdz  -  0  (4.27) 

-v2  -v2 

which  is  not  satisfied  by  the  term  mentioned.  However,  the  solvability  condition,  at  the 
next  order  requires  that  the  product  bQA1 1  is  always  zero,  which  alleviates  this  problem. 

Substitution  of  the  known  solutions  and  application  of  the  general  solvability 
condition  at  third  order  results  in  the  bifurcation  equations, 
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bo  Aii  “  0 


(4.28) 


‘  <a«2  '  0n2)T 
r  m  n 


*3b0  1A12  "  ° 


(4.29) 


which  along  with  the  normalizaiton  condition 


A1t2  +  A122  “  1 


(4.30) 


form  a  set  of  three  equations  for  the  three  unknowns  Aj  j,  Aj2*  bg.  The  term  a3  is 


given  by 


4  4  4  2 

5u)  3  (a  +8  )  3d). 

0  m  n  0 


- - 2_  + - ■  ..  —  ,  +  — L-  f|«  -  6  |  +  (a  +6)] 

3  32  „  4  ,,k  11  m  n1  m  n  J 


- -  "  ■■■-  [la  -  8  |  -  (a  +8)1 

_  /“  2  *•  m  n  m  m  J 

2/2  o.0 


2“02<2“02  *  V  2u02(2“02  "  V 


,  la  -  6  )  ( 3<i).4  +  8a  8  )  (a  +  8  )(3u 4  -  8a  8  ) 
l _  f  I  »  "I  0  -m  n  +  m  n  0  n  "  1 

_  L  .  2  i-  J 


32/2  [4«40  -/2|am-8j) 


[4u.^  -  /2(a_  +  8  )) 
u  m  n 


(4.31) 


It  has  not  been  proved  that  this  expression,  a3,  is  positive  definite  but  numerical 
evaluation,  with  m,n  ranging  from  1  to  10,  showed  this  to  be  the  case.  The  solutions 
to  (4.28)  -  (4.30)  are 


Case  It  b0  *  0,  A^  «  1,  A12  *  0 


/-(am2  -  Bn2)T 

Case  lit  b0  »  t  I - j -  .  An  -  0,  A12  -  1 

/  2*  3*0 


(4.33) 
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The  solution  is  (4.32)  is  of  little  interest.  It  merely  points  out  that  the  basic  solution 

bifurcates  from  the  primary  branch.  In  (4.33)  the  solution  for  the  secondary  bifurcation 

on  branch  1  is  given.  The  ±  sign  in  (4.33)  shows  that  the  bifurcation  takes  place 

in  both  the  upper  and  lower  half  planes  (±  £).  The  jump  to  the  *  0  solution  is 

often  referred  to  as  mode  jumping  because  the  solution  acquired  on  the  secondary  branch  is 

qualitatively  different  from  that  on  the  primary  branch.  The  radical  in  (4.33),  with  the 

conjecture  that  a3  >  0,  requires  that  (am  -  @m)T  <  0  for  secondary  bifurcation  to  occur 

on  branch  X 

m,n 

A  similar  analysis,  to  that  given  for  the  X  branch,  results  in  the  following 

m  t  n 

bifurcation  equations  for  the  X  branch, 

n  ,m 


(4.34) 

(4.35) 


which  along  with  (4.30)  provides  three  equations  for  the  three  unknowns.  For  convenience, 

define  e  to  be  the  point  of  secondary  bifurcation  on  the  X_  branch  and  e„  _  to 
n,n  m,n  n ,  m 

be  the  point  of  secondary  bifurcation  on  the  X  branch,  where 

n 

E  (U)  -  b  v  +  0(U2)  (4.36) 

m, n  n 

e  (u)  -  b  u  +  0(u2)  (4.37) 

n, m  n,m 

Equation  (4.33)  and  (4.34)  show  that  (retaining  only  the  positive  branch  for  brevity). 


The  secondary  bifurcation  phenomena  may  be  summarized  as  follows. 


V\,Vy?V‘V'.  wm 


I 


I 


i 


( 


k 


6 


I 

► 

► 

h 


P 


2  2V 

When  M  =  0  and  £  =  1  there  is  a  double  eigenvalue  A  «  (a  +8  )  /2 

0  m  n 


for  every 


pair  (m,n)  such  that  m  *  n.  As  £  is  varied  from  £  =  1  by  an  amount  £  =  1  +  Vjjni2, 
the  double  bifurcation  point  splits  into  two  primary  bifurcation  points. 


X  =  (a  2  +  £28  2)  1/2 
m,  n  m  n 


x  =  (8  2  +  £2a  2)  1/2 

n,m  n  m 


(4.40a) 

(4.40b) 


According  to  (4.36)  -  (4.39)  a  secondary  bifurcation  tkes  place  on  one  -  and  only  one  at  a 
time  -  branch.  Noting  that  T  =  sign  (£  -  1)  and  that  the  radical  in  (4.38)  and  (4.39) 
must  be  positive  for  secondary  bifurcation,  the  branch  on  which  the  secondary  bifurcation 
takes  place  is 


If 

£ 

> 

1 

and 

a 

> 

6 

then 

on 

X 

(4.41a) 

m 

n 

n,m 

If 

£ 

> 

1 

and 

a 

< 

6 

then 

on 

X 

(4.41b) 

m 

n 

m,n 

If 

£ 

< 

1 

and 

a 

> 

8 

then 

on 

X 

(4.41c) 

m 

n 

m,n 

If 

£ 

< 

1 

and 

a 

< 

8 

then 

on 

X 

(4.41d) 

m 

n 

n,m 

These  results  suggest  the  following.  As  £  departs  from  £  *  1,  the  split  primary 
bifurcation  points  move  away  from  the  double  point.  When  £  >  1  they  both  move  to  the 
right,  and  when  £  <  1  they  both  move  to  the  left.  However,  in  all  four  cases  (4.41),  the 
secondary  bifurcation  takes  place  on  the  branch  which  is  emitted,  after  splitting,  by  the 
largest,  in  magnitude,  of  the  two  bifurcation  points,  regardless  of  the  sign  of  T. 

Based  on  physical  and  numerical  grounds  the  size  of  the  £  neighborhood,  in  which 
secondary  bifurcation  may  take  place,  can  be  estimated.  Schwartz  and  Whitney  [10],  in 
their  analysis  of  two-dimensional  standing  waves,  in  a  fluid  of  infinite  depth,  estimated 
the  highest  wave  to  occur  at  e  ~  0.2  (the  scaling  here  is  different,  dividing  their 
result  by  n  provides  congruence  with  this  work).  This  two-dimensional  result  may  be  used 
as  an  estimate  of  the  relevant  three-dimensional  parameter.  Numerical  evaluation  of  the 
expression  for  bQ  gives  a  value  of  b^  ~  0.18  for  m  *  1,  n  »  2,  and  for  all 
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combinations  of  m,n  with  m,n  ranging  over  1  to  10,  the  value  of  bg  is  of  smaller 
magnitude.  This  suggests  an  upper  bound  for  p  of  p  ~  1.  Obviously  the  perturbation 
scheme  is  valid  only  as  p  ♦  0.  But  the  implication  is  that  the  neighborhood  around 
5  "  1  in  which  secondary  bifurcation  takes  place  may  be  substantial.  Therefore  the  range 
of  aspect  ratios  at  which  secondary  bifurcation  may  take  place,  at  some  value  of  the  wave 
amplitude,  is  large  enough  that  the  irregular  waves  produced  on  these  branches  should  occur 
quite  often  in  vessels  of  rectangular  cross-section.  In  the  next  section  the  solutions 
along  the  secondary  branches  will  be  found. 


\  wir*y  ■  jr  V  w?  w-9.  \*  77  t;  r;^ 
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5.  The  Solution  Along  the  Secondary  Branch 

After  the  splitting  process,  a  secondary  bifurcation  point  occurs  on  one  of  the  two 
resulting  primary  branches.  By  expanding  in  the  neighborhood  of  this  point,  an  asymptotic 
representation  of  the  solution  along  the  secondary  branch  may  be  found.  Assume  the 
parameters  are  such  that  the  secondary  bifurcation  point  occurs  on  the  Xm  n  branch.  A 
similar  expansion  can  be  performed  if  the  secondary  bifurcation  point  occurs  on  the  X^  m 
branch.  The  governing  equations  and  boundary  conditions  are  given  in  expressions  (2.1)  - 
(2.4).  The  variables  have  the  following  form, 


* 

»  e#  ip 

(5.  la) 

* 

n  =  eh  +  n 

(5.1b) 

* 

00  =  ui  +  ft 

(5.1c) 

where  4,  h. 

and  (1)  correspond  to  solutions  on 

the  primary  branch  and  $,  n. 

and  ft 

correspond  to 

solutions  on  the  secondary  branch. 

The  primary  branch  solutions 

were  found 

in  section  2. 

The  problem  for  1(1,  n,  and  ft  is 

given  by 

+  +  e2  =  0 

3x2  3y2  3z2 

(5.2a) 

at  x  -  ±  V2 

(5.2b) 

-  0  as  y  ♦  -  * 

(5.2c) 

*  0  at  z  =  ±  V2 

(5. 2d) 

and  on  y  =  eh  +  n 


*  01  M  +  rl a*  3,1  +  f2  3n  +  111  JLi  +  r2  111  JL1  l 

<M  +  ~  a?  +  3l  +  5  +  +  t  II  if  I 


+  9!  3n  +  2  li  an  =  .  n  3h  .  c2,3*  3h  *  2  3*  3h,  M 

3x  3x  3z  3z  3t  ’3x  3x  4  3z  3z>  3y 


(5.3a) 


and 


Ir 

► 
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« ♦  »>  M *  *  *  ■(££  * H ft  * i 


a*  e 

-E(o)  +  -  eh - -  V*«V$ 


The  free  surface  boundary  conditions  satisfied 
and  (5.3b)  after  expansion  of  these  conditions 
parameter,  v,  is  defined  as  a  measure  of  the 
point.  The  known  primary  branch  solutions  are 

e  -  b  w  +  0(u2) 
m,  n 

*  *  t0  +  P'P1  +  0(u2) 

h  “  hQ  +  uh1  +  0(y2) 

2  3 

u  =  id0  +  U(o2  +  0(u) 


by  ♦  and  h  are  subtracted  from  (5.3a) 
in  Taylor  series  about  y  =  0.  A  small 
distance  from  the  secondary  bifurcation 

(5.4a) 

(5.4b) 

(5.4c) 

<5.4d) 


with  bm>n  given  in  (4.38),  and  the  terms  in  (5.4c)  and  (5.4d)  are  given  in  (4.14a) 
through  (4.15b).  The  unknown  solutions  along  the  secondary  branch  are  postulated  to  have 
the  following  form 

♦  -  +  v2$2  +  .  .  .  (5.5a) 

n  -  vn1  +  v2n2  +  .  .  .  (5.5b) 

0  =  vf)1  +  v2S)2  +  .  .  .  (5.5c) 

Substituting  (5.4)  -  (5.5)  into  the  governing  equations  and  boundary  conditions,  (5.2)  and 

(5.3),  expanding  the  free  surface  boundary  conditions  in  a  Taylor  series  about  the 
equilibrium  solution,  and  equating  terms  proportional  to  like  powers  of  v  to  zero, 
results  in  a  sequence  of  boundary  value  problems  for  the  4^,  ,  and  However,  each 
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of  these  boundary  value  problems  is  also  a  function  of  p.  Therefore  the  functions 


(jf,  and  are  expanded  in  terms  proportional  to  powers  of  p. 


*i  =  H't'i,  +  V*  ^2  + 


ni  =  ur,i1  +  “  ^2  + 


n.  =  n.  n  +  pn. ,  +  p  +  .  .  . 

l  iO  ll  i2 


(5.6) 

(5.7) 

(5.8) 


The  substitution  of  these  expressions  results  in  an  additional  subdivision  of  boundary 
value  problems.  The  analysis,  although  straightforward,  is  lengthy  and  the  details  will 
not  be  presented. 


The  first  order  (in  v)  problem  results  in  Q1  *  0,  and 

2  3 

4'1  =  P^.,  +  ^  ^12  +  0(W  ) 

2  3, 

n1  =  pnn  +  p  n12  +  0(p  ) 


where 


t>. .  =  ii)„  ^ cos  8  x  cos  a  z  e 
r11  0  n  m 


n,«  =  cos  8  x  cos  a  z  sin  t 
ll  n  m 


“o  y 


cos  t 


(5.9) 

(5.10) 

(5.11) 

(5.12) 


where  the  parameters  are  as  previously  defined,  and 


4>. .  =  A  cos  (a  -  6  )x  cos  (a  -  8  )z  exp[/2|a  -  8  lyjsin  2t 
12  121  m  n  in  n  in  n 

+  A.-...cos  (<>  +  8  )x  costa  +  8  )*  exp[/2(a  +  8  )y]sin  2t 

122  m  n  m  n  in  n 

+  A  cos  (a  -  8  )x  costa  +  8  )*  exp(/2  w  2y]sin  2t 

123  m  n  m  n  0 

+  costa  +  8  )x  cos(a  -  8  )*  exp[/2  w  2y] sin  2t 

12  3  inn  inn  0 

and 


ni2  “ 

tB  1 2 1 
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C 1 2 1 

COS 
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8n)x 
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(am 

in 

6n)z 
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tB122 
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C122 
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(a 
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8  )x 
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cos 

(a 
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Bn)i 

+ 

tB123 
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C 1 2  3 

cos 
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cos 

(a 
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- 

B  )x 
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cos 

(a 
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tB123 
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C123 
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2t) 

cos 

(a 
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+ 

B  )x 
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cos 

(a 

m 

- 

6)  5 

m 

and  the  coefficients  in  $  and  n  _  are  given  in  Appendix  V. 


(5.13) 


(5.14) 


8u  [4u>  2  -  /l\ a  -  6  |  ]  8m  [4w  2  -  /2(a  +  B  )1 

0  0  1  m  n'  0  0  m  n 

+2i  "  ^1  “  °'  and 

26ny 

+22  “  A221  COB  26n*  *  sin  2t 
2a  y 

+  A  cos  2a  z  •  sin  2t 

222  m 

-  2“02y 

+  A__,  cos  2B  x  cos  2a  z  e  sin  2t  (5.19 

223  n  m 

and 

a 

n22  “  (B221  +  C221  CO“  2t)  008  2Cn* 

+  (B222  +  C222  «>•  2t)  cos  2aa* 

+  (B._,  +  cos  2t)  cos  2B  x  cos  2a  z  (5.20 

2 23  22 j  n  m 

The  coefficients  in  (5.19)  and  (5.20)  are  also  given  Appendix  V. 

The  result  (5.16)  provides  the  exprasslon  for  the  frequency  along  the  secondary 

branch.  In  the  limit  as  v  ♦  0  and  y  ♦  0  the  complete  expression  far  the  frequency  is 

*  2  2  2  3  3 

w  •  ♦  M  5^22  +  ,V  * 


(5.21 


Since  i8  proportional  to  quadratic  terms,  the  sign  of  determines  whether  the 

bifurcation  from  the  secondary  bifurcation  points  is  subcritical  or  supercritical.  It  has 

not  been  proved,  but  numerical  evaluation  of  the  expression  for  suggests  that  it  is 

negative  for  all  values  of  m  and  n.  Therefore,  as  the  solution  moves  from  the  primary 

to  the  secondary  branch  the  frequency  of  oscillation  decreases. 

An  example  of  this  secondary  bifurcation  phenomena  is  given  in  Figure  14.  In  this 

V, 

particular  example  m  m  2  and  n  «  1.  The  double  bifurcation  point,  Ag  2.65,  is 

2 

given  by  the  dashed  line  in  the  figure.  This  point  is  split  by  choosing  u  m  ^  .  With 
t  <  0  the  aspect  ratio  after  splitting  is  5  “  ^  •  Therefore  the  two  primary 
bifurcation  points  become  A.  .  ~  2.6  and  A  ~  2.4.  For  these  parameters  b_  ~  0.18. 

‘f  1  '»*  0 

* 

Therefore  the  secondary  bifurcation  takes  place  at  e  ~  0. 12,  on  the  A  branch,  and 

‘I  1 

from  this  point  the  subcritical  secondary  branches  are  emitted. 

The  nature  of  the  solutions  on  the  secondary  branch  can  be  illustrated  by  the  time 
series,  for  the  wave  height  in  the  left  front  corner  of  the  tank  (x  »  z  »  0),  given  in 
Figures  15a  and  15b.  The  result  corresponds  to  that  given  in  Figure  14.  In  Figure  15a 
the  wave  height  on  the  primary  branch  at  the  secondary  bifurcation  point  (e  ■  .12, 
v  »  0,  m  *»  2,  n  »  1 )  is  shown.  In  Figure  15b  the  solution  is  advanced  along  the 
secondary  branch  to  v  =  0.08. 

The  spatial  nature  of  the  solution  on  the  secondary  branches  will  be  similar  to  the 
result  found  for  the  mixed  modes  emitted  by  the  double  eigenvalues.  However  a  greater 
degree  of  complexity  may  be  expected  here  due  to  the  varying  amount  (through  v)  of  the 
"other"  mode  which  is  acquired  on  the  secondary  branch. 


A  bifurcation  diagram  with  the  occurrence  of  secondary 

bifurcation  for  m  «  2  and  n  »  1.  When  £  *  1  the 

pure  branch  emitted  by  the  double  eigenvalue  is  given 

by  the  dashed  line.  When  £  is  perturbed  to  £  *  7/9 

the  double  eigenvalue  splits  and  secondary  bifurcation 

occurs  on  the  A_  branch. 
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Figure  15a.  Times  series  for  the  wave  height  corresponding  to 

Figure  14.  The  parameters  cure  £  =  7/9,  m  «  2,  n  »  1, 
e  =  0.12,  and  vi  =  0.  This  is  the  solution  on  the 
primary  branch  near  the  secondary  bifurcation  point. 
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Figure  15b.  Time  series  for  the  wave  height  corresponding  to  Figure  14 
The  parameters  are  £  =  7/9,  m  »  2,  n  «=  1,  e  =  0.12,  and 
v  «  0.08.  This  is  a  solution  on  the  secondary  branch. 


6.  Remarks 


It  has  been  formally  shown  that  multiple  and  secondary  bifurcation  of  the  wave  fields 
occurs  in  partially  filled  vessels  of  rectangular  cross  section  when  the  depth  of  the 
fluid  is  infinite*  The  many  wave  solutions  depicted  in  this  paper  will  likely  occur  often 
in  nature.  Unlike  the  buckling  of  plates  and  shells,  for  example,  where  the  minimum 
eigenvalue  is  the  most  important,  in  the  "sloshing11  problem  the  band  of  frequencies,  at 
which  a  vessel  containing  fluid  might  be  excited,  is  relevant.  Consider,  for  example,  an 
ocean  liner  carrying  liquid  cargo  steaming  across  the  north  Atlantic.  The  ship,  and  hence 
the  cargo  tanks,  will  be  subjected  to  ocean  waves  containing  a  significant  band  of 
frequencies  and  amplitudes.  It  is  likely,  in  fact,  that  even  the  most  casual  observer  of 
waves  in  partially  filled  vessels  of  rectangular  cross-section  has  witnessed  wave  fields 
of  the  type  shown  herein. 

What  about  vessels  with  cross-sections  other  than  rectangular?  A  basin  with  a 
triangular  cross-section,  for  example,  will  have  a  multiple  (triple?)  eigenvalue  when  it 
is  equilateral.  Similarily,  one  can  imagine  that  a  star-shaped  cross-section  will  also 
produce  interesting  possibilities. 

A  basin  with  a  circular  cross-section,  such  as  a  coffee  cup,  is  however  different. 
There  is  no  relevant  parameter  such  as  an  aspect  ratio.  An  analogy  can  be  found  in  the 
work  of  Cheo  and  Reiss  [11]  on  the  secondary  buckling  of  circular  plates.  Their  analysis 
showed  that  a  different  type  of  secondary  bifurcation  occurs  in  the  buckling  of  plates 
with  circular  cross-section.  In  their  case,  the  primary  bifurcation  was  to  an 
axisymmetric  state,  from  which  a  secondary  bifurcation  to  an  asymmetric  state  was  found. 

In  the  standing  wave  problem  Mack  [12]  has  found  the  primary  bifurcation  points  to 
axisymmetric  states  for  standing  waves  in  a  circular  basin.  It  is  possible  that  a 
secondary  bifurcation  from  these  primary  states  to  asymmetric  states  My  take  place. 


The  nonlinear  boundary  conditions  at  the  free  surface,  y  ■  n  (x,*,t),  are 
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After  expansion  of  (1.1)  and  (1.2)  about  y  -  0,  and  substitution  of  (1.3),  the  following 
sequence  of  boundary  conditions  results 
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for  solvability 
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Appendix  IV 


In  this  appendix  the  coefficients  used  in  equations  (4.25)  and  (4.26)  for  the 


analysis  of  the  secondary  bifurcation,  are  defined.  They  are  not  to  be  confused  with 
those  in  Appendix  III. 
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In  this  appendix  the  coefficients  corresponding  to  the  solutions  along  the  secondary 


branches  are  given.  They  are  used  in  equations  (5.13),  (5.14),  (5.19),  and  (5.20). 
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